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Topological insulators are characterized by a nontrivial band topology driven by the spin-orbit 
coupling. To fully explore the fundamental science and application of topological insulators, material 
realization is indispensable. Here we predict, based on tight-binding modeling and first-principles 
calculations, that bilayers of perovskite-type transition-metal oxides grown along the [111] crystal- 
lographic axis are potential candidates for two-dimensional topological insulators. The topological 
band structure of these materials can be fine-tuned by changing dopant ions, substrates, and exter- 
nal gate voltages. We predict that LaAuOs bilayers have a topologically-nontrivial energy gap of 
about 0.15 eV, which is sufficiently large to realize the quantum spin-Hall effect at room tempera- 
ture. Intriguing phenomena, such as fractional quantum Hall effect, associated with the nearly- flat 
topologically-nontrivial bands found in e g systems are also discussed. 



Since the discovery^ 2 - of the quantum Hall effect 
(QHE), the quest for topologically ordered states of mat- 
ter has become a major subject of interest in condensed 
matter physics. Haldane first proposed 3 - that electrons 
hopping on a honeycomb lattice could realize the QHE 
in the absence of Landau levels, pointing out the pos- 
sibility of nontrivial topology in simple band insula- 
tors. Along this direction, recent efforts have culmi- 
nated in the theoretical prediction 4 ^— and subsequent ex- 
perimental realization^— of the so-called topological in- 
sulators (TIs) in materials with strong spin-orbit cou- 
pling (SOC). Many interesting phenomena, including gi- 
ant magneto-electric effects 1 ^ and the appearance of Ma- 
jorana fcrmions 1 ^, have been predicted. Once realized 
in real materials, these phenomena could lead to entirely 
new device paradigms for spintronics and quantum com- 
puting. 

However, so far, the material realization of TIs has 
been limited to narrow band-gap semiconductors based 
on Hg or Bi, in which the electronic properties are dom- 
inated by s and p orbitals. Here, we report our theo- 
retical investigation of topological insulating behavior in 
a completely different materials class — heterostructures 
of transition-metal oxides (TMOs) involving d electrons. 
Our motivation is two-fold. First, artificial heterostruc- 
tures of TMOs are becoming available owing to the re- 
cent development^ 3 . - — in the fields of oxide superlattices 
and oxide electronics 1 ^. In particular, layered structures 
of TMOs can be now prepared with atomic precision, 
thus offering a high degree of control over important ma- 
terial properties, such as lattice constant, carrier con- 
centration, spin-orbit coupling, and correlation strength. 
As we show below, these advantages can be readily ex- 
ploited in the design of TIs. Second, TMOs constitute 
a wide class of compounds that exhibit a variety of in- 
triguing properties and electronic states associated with 
the electron-electron interactions, encompassing super- 



conductivity, magnetism, fcrroelectricity and Mott insu- 
lators. Combined with the TI phase, TMO heterostruc- 
tures provide a very promising platform to explore vari- 
ous topological effects. 

Our main results are summarized below. We first 
demonstrate the design principle for realizing two- 
dimensional TIs in bilayers of perovskite-type TMOs 
grown along the [111] crystallographic axis by using phe- 
nomenological tight-binding modeling. Based on this de- 
sign principle and first-principles calculations, a number 
of candidate materials are identified. The topological 
band structure of these materials can be fine-tuned by 
changing dopant ions, substrates, and external gate volt- 
ages, which will enable also the control of the topologi- 
cal quantum phase transition. In particular, we predict 
that LaAuOs bilayer has a topologically nontrivial en- 
ergy gap about 0.15 eV, which is sufficiently large to re- 
alize the quantum spin-Hall effect at room temperature. 
When electron-electron interaction is included, our sys- 
tem with topologically nontrivial band structure could 
have far more interesting physics. Here, we demonstrate 
this possibility by focusing on the TMO bilayers of e g 
systems, which are characterized by nearly-flat topolog- 
ically nontrivial bands. We argue that when these 
bands arc partially filled, electron correlation could give 
rise to the quantum anomalous Hall effect and fractional 
quantum Hall effect. Our results may open new direc- 
tions focusing on topological phenomena in the rapidly 
growing field of oxide electronics. 
Results 

Design principle. To demonstrate the design prin- 
ciple for engineering TIs in TMO heterostructures, we 
consider perovskite-type TMOs as our prototype sys- 
tem. These compounds are very common and have the 
chemical formula ABO^, where O is oxygen and B is 
a transition-metal (TM) ion. The key idea is to start 
with a band structure that possesses "Dirac points" in 
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FIG. 1: Formation of the honeycomb lattice in a (111) bilayer in the cubic lattice, (a) Perovskite structure ABO3. 
(b) A (111) bilayer consisting of the top layer indicated by red circles and the bottom layer indicated by blue circles. The lattice 
constant is ao- The bilayer shown as solid lines in (b) forms the honeycomb lattice when projected on the [111] plane with 
the lattice constant a = ^/2/3ao (c). The real space coordinates are labeled by (x,y,z) in the original cubic lattice, while it is 
labeled by (X, Y) in the [111] plane, (d) Level structure of TM d orbital. In the cubic environment, d orbitals split into e 9 and 
t2 g manifolds. With the SOC, t2 g manifold further splits into two levels characterized by the effective total angular momentum 
j = 1/2 and 3/2. With the trigonal crystal field, t2 g manifold splits into two levels denoted by a\ g and e' g . With both the SOC 
and the trigonal field, t2 g manifold splits into three levels and e g manifold splits into two levels, i.e., all the degeneracies are 
lifted except the Kramers doublets, (e) ABO3 monolayer is grown on AO3 terminated AB' O3 substrate capped by AB' O3. 
The direction of crystal growth is indicated by an arrow. 



the Brillouin zone without the SOC, and then examine 
whether an energy gap can be opened at those points 
with the SOC turned on. If an energy gap does open, 
combined with proper filling the resulting state could 
be a TI. In an ideal perovskite structure, the TM ions 
sit on a simple cubic lattice, with the octahedral crys- 
talline field splitting the TM d orbitals into two-fold 
degenerate e g (c? 3z 2_ r 2, d x 2_ y 2) and three-fold degenerate 
t2g{d yz , d zx , d X y) levels, well separated by so-called IQDq 
on the order of 3 eV. Such a lattice geometry usually does 
not support Dirac points. Instead, we consider bilayers 
of the perovskite structure grown in the [111] direction. 
As shown in Fig. 1, the TM ions in the (lll)-bilayer are 
located on a honeycomb lattice consisting of two trig- 
onal sublattices on different layers. This lattice geom- 
etry has three consequences: Firstly, it is well known 
from the study of graphene that electrons hopping on 
a honeycomb lattice generally give rise to Dirac points 
in the band structure; secondly, A layer potential differ- 
ence can be easily created by applying a perpendicular 
electric field or by sandwiching the bilayer between two 
different substrates, which allows experimental control 
of the band topology; and, thirdly, the honeycomb lat- 
tice further reduces the symmetry of the crystalline field 
from octahedral (Oh) to trigonal (63,,), and introduces 
additional level splitting of the d-orbitals. The last point 
turns out to be crucial for realizing the topologically in- 
sulating phase. 

We first consider the t^g manifold, in which the on-site 
SOC is active. In our modeling, only nearest-neighbor 



hopping of d electrons between the TM sites via oxy- 
gen p orbital is included. Since we are interested in the 
band topology, which is robust against small perturba- 
tions as long as the band gap remains open, our model 
is justified and allows us to capture the essential ingredi- 
ents with minimal parameterization. The tight-binding 
Hamiltonian is given by 

H = -tJ2 T??'4r rf r'r' + lr ■ s r 
ty'tt' r 

+a 4^ + &4At , (i) 

where r and r label the lattice sites and the ti g orbitals, 
respectively. The first term is the hopping term rep- 
resented by a single amplitude t and the dimensionlcss 
structural factor T^i ■ The second term is the on-site 
SOC, which splits the ti g levels into a j = 1/2 doublet 
with energy A and a j = 3/2 quadruplet with energy 
—A/2. l r and s r are the angular momentum and spin 
operators. The third term is the trigonal crystalline field 
with which the ai g -e' g splitting is given by 3A/2. V in the 
last term is the layer potential difference, and £r = Tl 
when r is in the top or bottom layer. The explicit form 
of the Hamiltonian is presented in the Supplementary 
Method. 

The large number of orbitals (6 per TM site) involved 
in our model give rise to a very rich behavior of the topo- 
logical band structure in the parameter space. Depend- 
ing on the strength of the SOC, the system falls into two 
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FIG. 2: Dispersion relations of the (111) bilayer. (a) and (b) ti g model in the strong SOC limit. The SOC is fixed as 
X/t = 5 with A/t = 1 (red) and A = (green), (c) and (d) t 2g model in the weak SOC limit, A/t = 0.5 with X/t = 1.5 
(red) and A/t — 1.5 with X/t = (green), (e) and (f) e g model with X/t = 0.2 (red) and X/t = (green). Figures 
(a), (c), and (e) show the bulk dispersion relations. The dispersions in red correspond to the topologically nontrivial bands 
with the Z2 invariants shown for each band. Sum of Z2 in the occupied bands gives the Z2 topological invariant for the 
corresponding filling. For example, when the lowest 5 bands of the t 2g model are occupied by electrons in (a), Z 2 invariant 
becomes 1 + + + 1 + lmod2 = 1. The insets in (a) and (c) show the zoom-up near the K-point. Figures (b), (d), and (f) 
show the dispersion relations in finite-thick zigzag ribbons with the periodic boundary condition along the X direction and the 
openboundary condition along the Y direction. Parameters are the same as in the bulk dispersions. Edge modes supporting 
the spin current are indicated by red lines. For the t 2g model with the weak SOC, there appear 4 edge channels between the 
third and the fourth bands as shown as blue lines in consistent with the Z 2 number. 



different phases. In the strong SOC limit (|A/i| > 8/3 
when A = 0), bands originating from the j = 1/2 and 
j — 3/2 orbitals are completely separated. The trigonal 
crystal field then opens up an energy gap within each 
manifold and a nontrivial Zi topology can be realized. 
This is similar to the results reported on Iridium com- 
poundsi 17 ' 18 In the weak SOC limit, bands from j = 1/2 
and j = 3/2 orbitals become mixed away from the T 
point. Again, we find topologically nontrivial energy gaps 
can be opened by A. The Zi topological invariant is de- 
termined using two different methods. One can either 
evaluate it directly from the bulk band structure^ or 
count the number of edge states. The calculated band 
structure for both cases together with the Zi index is 



shown in Fig. [5] a-d. By inspection, we find that t\ g , 
t%i t~2g an d t% TMOs are all possible candidates for 
TIs in the strong SOC limit and t\ g and t\ g in the 
weak SOC. The dependence of the band topology on the 
layer potential difference V is rather interesting. While 
increasing V will eventually destroy the Z2 nontrivial 
phase, under moderate values of V ', the t| 9 system re- 
mains TI in the strong SOC limit, and the t\ g system in 
the weak SOC limit. 

Next we consider the e g manifold. It is well known 
that the SOC is quenched within the e g manifold so it 
seems that the resulting band topology should be trivial. 
However, similar to graphen o 20 ' 21 and some transition- 
metal ion s 22 ' 23 , the SOC can still take place through the 
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virtual excitation of electrons between e g and ti g levels. 
According to the second order perturbation theory, the 
effective SOC is given by 



zjes 

H SO ~ 



E 

r£e g 



(e\Hso\T)(T\H S o\e') 



(2) 



where e labels the e g orbitals. Hence the Hamiltonian 
can be written 



TT'ee' 



Ye 



£,vd Ve d Ye , 
(3) 

where the SOC magnitude is given by A = X 2 /A E , with 
A^- roughly being the energy difference between e g and 
tig levels. (The explicit form of the Hamiltonian is pre- 
sented in the Supplementary Method and Supplementary 
Note 1.) We find that Hsov opens up an energy gap at 
the r point and also the Dirac point located at K. From 
the inspection of the Zi topological invariant and count- 
ing the number of edge states, we found that e g , e 2 and 
e 3 systems become TIs (See Fig. [2] e and f ) . Here the 
trigonal crystalline field is also important — if all ti g lev- 
els are degenerate, even the second order SOC vanishes. 
A layer potential difference comparable to A closes a gap 
at the Dirac point turning the e 2 system into a trivial in- 
sulator. On the other hand, gaps at the T point are stable 
against this perturbation. Instead, these gaps close when 
the local potential difference between d 3z 2_ r 2 and d x 2_ y 2 
is comparable to A. Thus, the TI state and the Jahn- 
Teller effect^ compete in real materials with the e g or e 3 
configuration. 

Materials consideration. Having established that the 
TIs can be realized in (lll)-bilayer TMO for both ti g and 
e g configurations, we now turn to real materials. We aim 
to realize the integer fillings established above using TM 
B ions with the formal valence +3 or +4. For £? 3+ ( 4+ ), 
we choose La (Sr) for the A-site element in both the tar- 
get TMO and the insulating substrate AB'O^, and Al 
(Ti) for the B'-site element in the insulating substrate. 
Controlling the strain effects and the layer potential dif- 
ference is possible by replacing A and/or B' with their 
isovalent elements. It is well known that some of TMOs 
are insulating due to strong correlations^. Therefore, if 
the corresponding bulk system is heavily insulating, bi- 
layering may not be useful. Even if the corresponding 
bulk system is metallic, the low dimensionality in (III) 
bilayers may drive the system into a Mott insulator- 6 . 
Further, the correlation effects are expected to reduce 
the effective band width and increase the splitting be- 
tween occupied levels and unoccupied levels. While this 
effect does not change the band topology in e g electron 
systems, this could influence the topology in ti g systems 
by modifying the crystal field splitting between a\ g and 
e' g levels. In addition, in a system with an integer-number 
of electrons per site, local moments could be induced by 
the correlation effects resulting in the magnetic ordering. 
If the symmetry breaking by the magnetic ordering is 
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strong, the system could become a trivial insulator. We 
do not consider such complexities by focusing on rather 
itinerant 4c? and 5d electrons of TM ion, yet ti g electron 
systems are more susceptible for magnetic orderings than 
e g electron systems because of the smaller hopping inten- 
sity. These considerations somewhat limit the choice of 
TM and substrate material. Our candidate materials for 
TIs are, therefore, LaRe 3+ 03 as a t\ electron system, 
LaRu 3+ 3 , LaOs 3+ 3 , SrRh 4 +0 3 and SrIr 4 +0 3 
systems, and LaAg 3+ 3 and LaAu 3+ 3 as e 2 
systems. Most of these materials have been synthesized 
and their references are summarized in Table HI LaRe0 3 , 
LaOs0 3 and LaAg0 3 have yet to be synthesized. Accord- 
ing to Ref. HH, LaAu0 3 has CaF2 structure rather than 
the perovskite. We expect this material shapes the per- 
ovskite structure by, for example, high-pressure synthesis 
and grown on a substrate with the perovskite structure. 
If properly synthesized, perovskite LaAu0 3 is expected 
to be metallic as LaAg0 3 predicted by the density func- 
tional theory calculation^ 3 -. We also note that growing 
thin films of perovskite TMOs along the [111] direction 
has already starte d 36 ' 37 . While t 2 g systems are also can- 
didates for TIs, the TI state is hard to realize because 
of the band overlap. For e l g and e 3 systems, additional 
effects such as longer range transfer and the Jahn- Teller 
effect can easily modify the dispersion relations. 

We first performed the density functional theory 
(DFT) calculations for the bilayers of ti g systems 
LaRe0 3 , LaRu0 3 , LaOs0 3 , SrRh0 3 , Srlr0 3 . (Details 
are presented in the Method section.) Their dispersion 
relations are shown in Fig. [3] a-d. We notice the re- 
markable agreement between the DFT results and the 
tight-binding result, Fig. [5] c, especially for LaRe0 3 
and LaOs0 3 . For SrRh0 3 , LaRe0 3 , and LaRu0 3 (not 
shown), the Fermi level crosses several bands. Thus, 
these systems are classified as topological metals rather 
than TIs. In LaOs0 3 and Srlr0 3 , the Fermi level is lo- 
cated inside the gap. Therefore, from the analogy to 
the tight-binding model, (111) bilayers of LaOs0 3 and 
Srlr0 3 are TIs. From our DFT calculations, it is noted 
that the material dependence of the dispersion relations 
is rather large for ti g systems. This is because a large 
number of band parameters arc involved in the band 
structure including the local crystalline field. 

We now move to e 2 g electron systems, LaAg0 3 and 
LaAu0 3 . Our DFT results for these systems are shown 
in Fig.[3]e and f. As in the ti g case, the DFT reproduces 
the tight-binding result fairly well. In both undoped sys- 
tems, the Fermi level is inside the gap at the K point, 
and from the analogy to the tight-binding result, these 
systems are TIs. The gap amplitude is found to be about 
150 meV for LaAu0 3 and 40 meV for LaAg0 3 , so these 
systems should remain TIs at room temperature. The 
band gap and topological property can be controlled by 
breaking the symmetry between top and bottom layers. 
As shown in Fig.Og, the asymmetric bilaycr with LaSc0 3 
has larger band gap ~ 300 meV, and from the inspection 
of the symmetry of the wave function at the K point, 
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TABLE I: List of candidate materials 



Configuration 



Bulk 



Superlattice 



LaReC>3 


t 4 


LaRuOa 




SrRhOa 


r' 


SrlrOa 


r' 


LaOsOs 




LaAg0 3 


12 

e s 


LaAu0 3 


4 



metallic Ref. 27 
metallic Ref. .28 
metallic Refs. .30 & .31. 

metallic (band calc.) Ref. 33 
Refs. 34 & 35 



Ref. 29 
metallic Ref. 32 




FIG. 3: Density functional theory results of the dispersion relations of the (111) bilayer of TMOs. Symmetric 
bilayers: (a) LaRe03, (b) LaOsOs, (c) SrRh03, (d) SrlrOs, (e) LaAgOs, and (f) LaAu03. Bilayers shown in (a), (b), (e), 
and (f) are grown between LaA103, while those in (c) and (d) are grown between SrTi03. Asymmetric bilayers of LaAuC>3 
grown between LaAlC>3 and LaSc03 (g), and between LaA103 and YAIO3 (h). The Fermi level is taken to be of the vertical 
axis. Bilayers shown in (b), (d), (e), (f), and (h) are TIs with the band gap indicated, (g) is a trivial insulator, and others are 
topological metals. 
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this bilayer is a trivial insulator. While the asymmetric 
bilayer with YAIO3 has smaller band gap ~ 50 meV and 
remains to be a TI (Fig. [3]h). Such a small band-gap TI 
is especially useful to control the topological property by 
using the gate voltage. Similar to LaAgC>3 and LaAuC>3, 
we also performed the DFT calculations for a 3d system 
LaCuC>3. We found this system develops an instability 
towards magnetic ordering because the itinerancy is re- 
duced compared with 4d and 5d systems. 
Nearly-flat Z 2 topologically nontrivial bands. One 
of the appealing aspect of realizing nontrivial band topol- 
ogy in TMOs is the rich possibilities of novel phenomena 
that could emerge when electron correlation is consid- 
ered. Here we demonstrate one of the possibilities in the 
e g systems with nearly-flat Z2 topologically nontrivial 
bands (see Fig. [5]c). Obviously, when the chemical po- 
tential is tuned into the nearly-flat Z2 bands, we are in a 
novel regime of quantum many-body physics. From the 
inspection of dispersion relations shown in Fig. [3]c and f, 
the upper fiat band appears to be more stable than the 
lower flat band. Tuning the chemical potential in the up- 
per flat band corresponds to removing electrons from d 10 
systems or adding electrons to d 9 systems. When such 
a situation is realized, kinetic energy is suppressed and 
physics is controlled mainly by interaction, whose energy 
scale is typically 1-2 eV, much larger than the width of 
flat bands W ~ 0.2 eV. What would be the ground state 
of the system? Here we propose several natural candidate 
states. One very likely consequence of the short-range re- 
pulsion U, when the doping of the nearly-flat band is not 
too small, is to drive the system ferromagnetic because 
of Stoner's criteria U /W ^> 1. While when doping is too 
small, Wigner crystal phase should naturally occur. In 
the following we assume that spontaneous ferromagnetic 
ordering occurs. The magnetic order should be viewed 
as a breaking of a discrete symmetry, due to the SOC, 
instead of the breaking of a continuous spin rotation sym- 
metry which is not realized in a two-dimensional system 
at any finite temperature. Our discussion is not limited 
to spontaneous ordering, because magnetism can be also 
introduced externally by using a magnetic insulator as a 
substrate. 

We model the effect of ferromagnetic ordering by 
adding a Zeeman spin splitting to the Hamiltonian: 

Hzeeman = — B ^^(d r ^.dret — 'Aei^^i) '■ (^) 
Ye 

Resulting dispersion relations for e g system with two 
characteristic Zeeman fields are presented in Fig. [5] From 
the inspection of the Chcrn number, the quantum anoma- 
lous Hall (QAH) insulating stated is realized in e°' 5 and 
e^- 5 systems, and also e x g systems when the Zeeman field 
is large as in perovskite manganites. More interesting 
physics occurs when the nearly-flat "Chern" band is par- 
tially filled (See Figs. U e and 13} . In this pointed 
out recently, fractional quantum Hall (FQH) states are 
likely to be realized^— . To elaborate this possibility, 
we perform exact diagonalization calculation after pro- 



jecting on-site repulsion and nearest-neighbor repulsion 
into the 1/3 filled highest Chern band (i.e. at eg' 5+1 ^ 6 ). 
Indeed, signatures of a v = 1/3 FQH state are observed 
(the detail of the exact diagonalization and the numerical 
result are presented in the Supplementary Note 2). We 
find that, when the nearest-neighbor repulsive interac- 
tion is larger than the width of the flat band, the ground 
state degeneracy is three-fold on a torus, and the Chcrn 
number (excluding the integer Chern number from the 
filled bands) of the ground state wave function is ~ 1/3 
up to finite size correction. 

The aforementioned FQH effects and QAH effects are 
both high-temperature effects and fundamentally dif- 
ferent from the quantum Hall states realized in GaAs 
two-dimensional electron gas in a magnetic field, where 
quasi-particle energy gap is controlled by the long-range 
Coulomb repulsion e 2 /(el B ) (e is the dielectric constant 
and Ib is the magnetic length), typically around a few 
Kelvin. In the present systems, quasi-particle energy gap 
is determined by the short-range repulsion ~ 1-2 eV. 
This indicates that room-temperature FQH effects may 
be realized. FQH states, in particular the non-Abclian 
states, have been shown to be very useful as building 
blocks of a quantum computer—^. A high tempera- 
ture non-Abelian quantum Hall states in the TMO het- 
erostructures, for example, at v = 1/2 filling where nat- 
ural candidate states are in the same universality class of 
Pfaffian states^ or anti-Pfaffian state a 48 i 49 , if realized ex- 
perimentally, would have strong impacts on both funda- 
mental physics and its applications, including the efforts 
of realizing topological quantum computation. 
Discussion 

Before closing, wc make few remarks on TIs in the TMO 
bilayers. The direct confirmation of the TI state is pos- 
sible by measuring the conductance. As in Refs. [H and 
8', the conductance should be quantized as a — 2e 2 /h 
per (111) bilayer in the two-terminal measurement. The 
conductance can be controlled by using the gate voltage. 
In our DFT calculations, the (111) bilayers are repeated 
along the [111] axis. With the non-zero inter-bilayer cou- 
pling, the helical edge channels on the surface of the sam- 
ple will turn to the two Dirac fermions at fc[m] = and 
1/2 in the unit of the reciprocal lattice vector along the 
[111] direction. Thus, strictly speaking, the TI is classi- 
fied as a "weak" TI and the backward scattering between 
the two Dirac fermions, if it exists, causes the localiza- 
tion. In order to avoid this, one needs to make the single 
(111) bilayer or keep the neighboring (111) bilayers far 
apart so that the inter-bilayer coupling becomes exponen- 
tially small. When fabricating the (111) bilayer, a small 
number of defects would have minor effects because the 
edge modes surrounding them are disconnected. How- 
ever, as the defect density increases, two surfaces are 
eventually connected through the edge modes belonging 
to the islands. As a result, the backward scattering takes 
place. Another source of the localization is islands of 
thicker (111) TMO layers. The suppressed trigonal field 
inside such thick islands is expected to create the gaplcss 
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FIG. 4: Bulk dispersion relations of the (111) bilayer of e g model with the Zeeman splitting, (a) Small Zeeman 
splitting with B = 0.3t, and (b) large Zeeman splitting with B = 2t. Here we used A = 0.5t. Bands with the majority 
(minority) spin component are indicated as red (blue) lines. Band dependent Chern number is also indicated. 



bulk modes. In further thicker islands of (111) TMO lay- 
ers, metallic state would be realized inside the sample. In 
addition to (111) bilayers, we have studied model (111) 
trilayers and found that the TI states are robust without 
the even-odd oscillation between TI and trivial insulator, 
as predicted for the bismuth thin films^. Structurally, a 
(111) trilayer forms a so-called dice lattice, which could 
also bring about interesting quantum effects character- 
ized by the Chern number C = ±2^1. Of course, if the 
layer structure is too thick then the bulk cubic symme- 
try is restored and the system is no longer a TI. So far, 
we did not mention the correlation effects and compet- 
ing ground states except for the previous section. For 
limiting cases, we have performed unrestricted Hartree- 
Fock calculations for multi-orbital Hubbard models de- 
fined on (111) bilayers. We found that the TI states 
are rather robust for e g systems and become unstable 
against antiferromagnetic insulating states when the in- 
teraction strength is comparable to the full bandwidth 
as in the two-dimensional Hubbard model on the honey- 
comb lattice^. For e g or e g systems, the QAH insulating 
states could be generated dynamically by correlation ef- 
fects without the SO C 53 ' 54 , yet trivial insulating states 
due to the Jahn- Teller effect would also be stabilized de- 
pending on the relative balance between the Coulomb 
interaction and the Jahn- Teller coupling. In this paper, 
we focused on the perovskite-type TMOs. Thus, our de- 
sign principle for the TI state works only for the [111] 
plane because other planes such as [001] and [110] do not 
support a honeycomb lattice. However, this approach is 
not limited to the perovskite systems. For example, the 
[0001] plane of corundum AI2O3, i.e., sapphire, involves 
a honeycomb lattice formed by Al atoms. Such a sys- 
tem could also be utilized as the substrate material to 
artificially create the TI state. 
Methods 

Tight-binding models in the real space. First, we 
consider a general multiband tight-binding (TB) model 
on a cubic lattice given by 



Hband 



(rr')a MM' 



•}. (5) 



where r labels the transition-metal sites, a spin, and fi 
orbitals. is a transfer matrix which depends on the 
pair of orbitals but not on the spin; its detail will be 
presented shortly 

For t2 g electron systems, the trigonal crystal field di- 
rectly couples with the local ti g level. In addition, the 
angular momentum is not quenched, and therefore the 
spin-orbit coupling (SOC) is active. Including these 
two effects, a TB model for t^g systems is written as 
H t 2g = H band + H S o + H tri with H S o and H tri given by 
the second and the third terms of Eq. (1), respectively. 
The explicit form of Hso f° r the £2g-al°ne model is given 
by 



Hso = A E lr ' Sr = ^E E i£TT'T"d\. Ta a T cnyl d VT , <y , 

Y Yaa' tt't" 

(6) 

with the use of the following convention for the orbital 
index: |o) = \d yz ), \b) = \d zx ), and |c) = \d xy ). a T with 
t = a, b, c is the Pauli matrix, and e TT > T " is the Levi- 
Civita antisymmetric tensor. 

The dependence of transfer matrices on the orbital and 
direction is given by the Slater-Koster formula^ as fol- 
lows: 



_ .66 

"r,r±y(z) r,r±z(x) T,r±x(y) 



t" 



= f 



hi 



r,r±x r,r±y r,r±z 
for the nearest-neighbor (NN) hopping and 



iflfl j.66 tec 1 

i\r±y±z — r,r±z±x — r,r±y±z ~ 



(7) 
(8) 

(9) 
(10) 



j.ab j.oc j.ca ± 

r,r±(x+y) " r,r±(y+z) ~ r,r±(z+x) ~ 

t ab - - — t bc - - - t ca - - = -t 1 (11) 
r,r±(x-y) [ T,r±(y-z) t r,r±(z-x) b7T > y ±x ) 

for the second-neighbor (SN) hopping. Here, x, y, and 
z are the unit vector along the x-, y-, and z-direction, 
respectively. Although it is via weak it hybridization t^ d 
between a transition-metal ion and an oxygen ion, the NN 
hopping t v oc (tp d ) 2 /A pd is the largest parameter in this 
model, thus, taken as the unit of energy t. A pc i is the level 
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difference between TM d orbitais and oxygen p orbitals. 
The ratio between and t v is the dimensionless param- 
eter T^yi = ^rr'/^Tr- ts> is a l so the NN hopping, but it is 
via weak direct overlap and, therefore, is expected to be 
small. t a ii and t v i are the SN hoppings due to the higher- 
order processes involving the transfer between two oxy- 
gen ions as t a » ^' °c {tpd) 2t pp I '^>d- Since involves 
relatively strong (weak) a(jr) hybridization between two 

oxygen ions tpp , we expect |t CT «| > \t w '\. Typical trans- 
fer intensities are shown in Supplementary Figure SI (a). 

For e g electron systems, linear coupling with the trigo- 
nal crystal field is absent. Therefore, the C3 lattice sym- 
metry of a (111) bilayer does not influence the on-site 
e g level. On the other hand, e g degeneracy can be lifted 
by the distortion of an 06 cage surrounding a transition- 
metal ion, i.e., the Jahn- Teller effect. Focusing on the 
metallic regime, we neglect this effect. The angular mo- 
mentum is quenched unless the coupling between e g and 



Here, s(= a,f3) labels the e g orbitals as \a) = \d 3z 2_ r 2) 
and |/3) = \d x 2_ y 2). For e g electron systems, the NN 
hopping t a oc (tp d ) 2 / Apd is via strong a hybridization 
between a transition-metal ion and an oxygen ion t° d 
and, therefore, largest. This hopping integral is taken as 
the unit of energy t. Again, the dimensionless parameter 
T^yi is defined by the ratio between t^ E r , and t a as T^, = 
^fr'Ac- The NN hopping ts is due mainly to the direct 
overlap between two TM ions and, therefore, expected 
to be small \t$\ -C \t a \ as tgi in the t2 g orbital model. 
t a * is the SN hopping due to the higher-order processes 
involving the transfer between two oxygen ions as t a > cx 
(tp d ) 2 tp p / Ap d . Typical transfer intensities are shown in 
Supplementary Figure SI (b). 

Tight-binding models on the (111) bilayer. By 

constraining the atomic coordinate r within the (111) 
bilayer R = (X, Y), it is straightforward to derive the 
TB Hamiltonian function of two-dimensional mo- 



t 2g orbitals are considered. We also neglect this effect mentum k = (k x ,k y ). We use the convention in which 



at the moment but reconsider it later. Thus, for the e g 
alone model, H e = Hband- The dependence of transfer 
matrices on the orbital and direction is again given by the 
Slater-Koster formula^. For the NN hopping, we have 



r,r±z 



r,r±x 



r,r±y 



j.aa j.' 

r.r±z r,r±z 

,0a 



,ct0 

r,r. 





X 0" 




.0 t s . 



t a + 3t 4 -V3(t a - t 5 ) 
-V3(t a -t s ) 3t a + t s 

t a + 3ts V3(t a — ts) 

V3(ta — ts) iter + ts 



(12) 
(13) 
(14) 



the projection of the nearest-neighbor bond into the 
[111] plane a is taken as the unit of the length scale. 
This is a factor \/2/3 smaller than the lattice constant 
of the cubic perovskite, and the size of the new unit 
cell is 3\/35 2 /2. Taking the primitive lattice vectors 
as ai = (\/3d,0) and &2 = (v / 3a/2, 3a/2), the first 
Brillouin zone is a hexagon with 6 corners located at 
k= (±47r/3\^a,0),(±27r/3\/3o,±27r/3a). 



and for the SN hopping 



t 



r,r±x±y 



1 

2*" 



t ee ~ = —-t 
r,r±z±x 2 



t ee . — + 

r,r±z±y 2 



1 
-3 j 

2 -V3 
-a/3 

2 x/3 
V3 



(15) 

(16) 

(17) 
I 



It is straightforward to derive the Hamiltonian matrix 
in the (111) bilayer. For the t2 g orbital model, we obtain 



band 



Htri — 



k 



lak 


1 


d lk 




ick 




2ak 




4k 




-^2ck- 





V _1_ B 


^ck 
V 1 - 
~T+ £ bk 


%k 


£ ak 






%k 


lak 

V _L B 




£ bk 




7 h k 


7„.k 






£ rk 








\+hk 


T c k 


%k 




£ lk 




7 c k 


V 1 B 
-+ £ bk 


Tak 






F* 
£ ck 


%k 


7 a k 


T+ £ ck. 





d lak 




d uk 




d ick 




d 2ak 




d 2bk 




- d 2ck - 



with 



~ak 



.{l 



-Ml 



V3r 



ts>e 



hk 



19) 



-t n {l 



9 



£ ck 

7ak 

%k 

%k 

£ ak 

£ bk 
£ ck 



-t s > 



1 



- 2^ COS^ — fc^ 

/ a/3, 
2^< cos(^ — — fc^ 

2i 7r < COS x/Sfcj;, 

——k x + 



3fc 



-A + 2t OT / cos 



-2t„ 
-2t c 

-2tr. 



h 

2 v 



h 
2 y 



COS K T 



2 

cos \/3fc.T 



2*" 



(21) 

(22) 

(23) 
(24) 

(25) 

(26) 
(27) 



Here, the spin indices are suppressed for simplicity. 
is the sublattice-dependent potential which breaks the 
symmetry between the top (labeled 1) and bottom (la- 
beled 2) layers. 



For the e„ orbital model, we obtain 



band 



E 

k 



lak 


T 


r v , ~ 

~2 + £ a k 


£ aflk 


u 




£ apk 




2ak 




ok 


afik 


d f u 






F* 



rfk 



£ ak 
£ a/3k 
^+ £ a k 
£ a/3k 



£ a pk 

£ fik 
£ a pk 

r + £ pk 





d l a k 




d l/3k 




d 2ak 




dzpk- 



(28) 



with 



£ a k = -t* - \% + 3t s) cos (jY k * ^ 



V3, 



£ [3k = -** - 2^ ta + t $) C0S y : Y k * 

v3 ,. . \ ■ . ( Vo , 
£ a pk = ~ — (U-U)isml — k x )e 



(29) 
(30) 
(31) 



£ a k = ta>^cos(-^-k x \ cosQfcy^ -cos\/3fc x |, (32) 

£ sk = 3*o-' cos V^fc^, (33) 
_ /a/3 \ /3 \ 

(34) 



£ a pk = 2 a/3V sin( sm(^k v 



DFT calculations. Density functional calculations 
were carried out using the projector augmented wave 
(PAW) method^ with the generalized gradient approxi- 



mation (GGA) in the parametrization of Pcrdcw, Burke, 
and Enzerhof (PBE)££ for exchange-correlation as im- 
plemented in the Vienna Ab Initio Simulation Package 
(VASP) 5? . The default plane- wave energy cutoff for 
O, 400.0 eV, was consistently used in all the calcula- 
tions. The optimized crystal parameter is 3.81 A for 
bulk LaAlCb and 3.95 A for bulk SrTi0 3 . These values 
are in consistent with the exp erimental values of 3.79 A 
(Ref HI) and 3.91 A(Ref. M). The transition-metal bi- 
layer structures were simulated by a supercell consisting 
of 12 AOs and 12 B layers along the [111] direction with 
(A, .B)=(La,Al) or (Sr,Ti) with two adjacent B layers re- 
placed by transition-metal ions. In the (111) plane, the 
supercell contains a 1 x 1 unit cell. A 6 x 6 x 1 special 
k-point mesh including the T point (0,0,0) was used for 
integration over the Brillouin zone. Optimized atomic 
structures were achieved when forces on all the atoms 
were less than 0.01 eV/A. 
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Supplementary Figure SI. Orbital dependent transfer integral on a (111) bi- 
layer for (a) t^g model and (b) e g model. Wave functions and lattice sites on the 
top layer arc indicated by red and those on the bottom layer arc indicated by blue. 
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Supplementary Figure S2. Energy spectrum of the effective model Eq. f|S2.2[) 

defined on a 4 x 6 cluster with 8 electrons. Nearly degenerate GSM and the 
large gap between GSM and excited states indicate the formation of the fractional 
quantum Hall states. 
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Supplementary Notes 

Supplementary Note 1. Effective spin-orbit coupling for e g system in the (111) bilayer 

For an e g model in the cubic symmetry, the SOC is absent because the angular momentum is quenched. However 
in (111) bilayers, the local symmetry is lowered to C3 and, therefore, the virtual coupling with ti g orbital can induce 
the effective SOC between the two e g orbitals. 

To see this effect, we consider the a\ g orbital given by the linear combination of three tig orbitals. This state is 
expected to be highest in energy among three orthogonal orbitals originating from ti g triplets. Two e g orbitals are 
expressed as \a) = \d 3z 2_ r 2) = \m = 0) and \/3) = \d x 2_ y 2) = — ^ { 1 2) + | — 2)}, where the azimuthal quantum number 
I = 2 is suppressed for simplicity. Here, the spin quantization axis is taken to be z in the original cubic lattice. The 
wave function for a\ g is expressed as 



= i={-(l - i)\l) + (1 + i)\ - 1) - i(|2) - I - 2))}. 



(Sl.l) 



Using the atomic SOC, Hso = Al • s = \{l z s z + i(7 + s_ + /_s + )}, as a perturbation, and the level separation 
between e g and ai g , we obtain the effective SOC for the e g -alone model as 



The diagonal term X/2Iq is absorbed in the chemical potential. The effective SOC constant is given by A = X 2 /Ae- 
Therefore, when the level separation between e g and a\ g orbitals is small, the effective SOC A could become large. 
Other contributions including the excitation to e' g orbitals are expected to normalize A. In particular, when a\ g 
and e' g orbitals are degenerate, i.e., A = 0, off-diagonal components of Eq. (|S1.2[) disappear. By changing the spin 
quantization axis along the [111] direction, the effective SOC is conveniently written as 

E^or = -Ei^C^v, (S1.3) 
r r 

where t v is the Pauli matrix acting in the orbital space. 

Supplementary Note 2. v — 1/3 Fractional Quantum Hall state on e g flat band 

The tight-binding model that we are studying is given by Eqs. (28) and (|S1.2|) . We choose A = 0.5i, B = 0.3i, and 
the 8-bands are shown in Fig. 4 a. 

Now we focus on the highest band - the 8th band (it differs from the lowest band by a particle hole transformation) , 
and consider the electron interactions Hi. For Hi, we use the following interaction: 

Hi = U n Te ^n Tei + U' ^ n r e "r e ' + V ^ n T n T '. (S2.1) 
re re>e' <rr') 



-'rat 
u rai 



1 + i 



1 - i 





drat 




drpt 




drai 




drpi 



(S1.2) 



Here r labels site, e labels e g orbitals, and a =f, I labels spins. nr ecr = d Yea dv e cr is the electron density for orbital-e 
and spin-CT, nr £ = nr £ a is the electron density for orbital-e, and = Y] £ n-r £ is the total electron density at 
site-r. U is on-site intraorbital repulsion, U' is on-site intcrorbital repulsion, and V is nearest-neighbor repulsion. 

Now we 1/3-fill the 8th band, and try to find out the ground state. We project the interaction Hi into the 8th 
band, and study the effective Hamiltonian in the partially filled band (the 8th band) . This treatment does not include 
the band- mixing due to the interaction so it is not exact. But it is a reliable treatment when interaction U, [/', 
V are weak compared with the band energy spacing. In addition, band-mixing can be included in future study by 
perturbation treatment. 
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H eff = Y^Es(k)^ k ^ + — ]T "(ki.ka.kaK^^^+k,-^. ( S2 - 2 ) 
k uc kik 2 k 3 

where Eg(k) is the kinetic energy of the 8th band, and the interaction u is nothing but Hj projected into the 8th 
band. N uc = N x ■ N y is the total number of unit cells. 1/N UC is the correct normalization factor. 

H e jf is exact diagonalized for 4 x 6 (N x x N y , these are the number of unit cells along a! and a 2 directions on 
honeycomb lattice) unit-cell system with the periodic boundary condition, with N e = 8 electrons. Because H e ff 
respects the total momentum, we can diagonalize H e ff with in each center-of-mass (COM) momentum sector. In 
Supplementary Figure S2, we show the ground state and the first excited state energies as a function of the center- 
of-mass momentum. We choose U = V = t,V = 0.5t. The momentum k x ,k y are shown as integers. For example, 
(k x , k y ) = (2, 3) really means (k x , k y ) — (2 • 2w/N x , 3 • 2tt/N v ). 

A three-fold degenerate ground state manifold (GSM) is observed, which is separated with the other states by a clear 
energy gap ~ 0.1. Note that they are not exactly degenerate on a finite system. But, their energy difference should fall 
off exponentially as the system size increases. These three ground states are at momentum (0,0), (0,2), (0,4) which 
is expected. The reason is that the different ground states can be viewed as a result of twisted boundary condition 
— > 2n. If we twist the boundary condition along the y direction — > 2tt, the momentum of each electron is shifted: 
k y — >• ky + =^-, So for 8 electrons, the center of mass momentum shifts k y k y + 8^ = k y + This twist will 
drive ground state 1 with center of mass momentum (0,0) to ground state 2 with COM k = (0, 2). And it also drives 
ground state 2 into ground state 3. 

To confirm that this 3-fold degenerate ground state is really a fractional quantum Hall (FQH) state instead of 
states such as CDW, we computed the Chern number by twist boundary condition. The details of the method are 
described in Ref. 61. Here we discrctize the boundary phase unit cell into a 10 x 10 and 20 x 20 meshes. And the 
Chern numbers for the three ground states are found to be independent of which mesh to use up to the fourth digit: 
C\ = 0.3344, C2 = 0.3311, C3 = 0.3344. These values slightly deviate from C = 1/3 in the thermodynamic limit, 
which is expected for a small system. The sum of the three Chern numbers is found to be exactly 1. This explicitly 
shows that we are in the v = 1/3 FQH phase. 
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